Creating Ground State Molecules with Optical Feshbach Resonances in Tight Traps 
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We propose to create ultracold ground state molecules in an atomic Bose-Einstein condensate by 
adiabatic crossing of an optical Feshbach resonance. We envision a scheme where the laser intensity 
and possibly also frequency are linearly ramped over the resonance. Our calculations for *^Rb 
show that for sufficiently tight traps it is possible to avoid spontaneous emission while retaining 
adiabaticity, and conversion efficiencies of up to 50% can be expected. 

PACS numbers: 32.80.Qk, 33.80.-b, 34.50.Rk 



The formation of ultracold molecules and the creation 
of molecular Bose-Einstein condensates (BEC) open 
the way to study new collective phenomena and a new, ul- 
tracold chemistry 0, Q . Since no direct cooling method 
for molecules can reach the transition temperature for 
BEC, the formation of molecules from ultracold atoms 
has been a focus of recent research. Molecules are cre- 
ated by applying an external field, either magnetic 
or optical pj, to two colliding atoms. This process is 
described in terms of a Feshbach resonance (FR) 6] 
where the collision energy of the two atoms coincides 
with the energy of a bound molecular level. Magnetic FR 
have been particularly successful in creating alkali dimer 
molecules 0,0], even heteronuclear 0. In contrast, op- 
tical FR involve electronically excited potentials, where 
spontaneous emission may lead to loss of coherence 0,0 . 
Apart from this obstacle, optical FR have the advan- 
tage that optical transitions are almost always available, 
whereas magnetic FR require the presence of a hyperfine 
manifold of the atom and may occur at magnetic field 
strengths which are difficult to obtain in experiments. 
Furthermore, optical FR offer more flexibility since two 
parameters (laser intensity and frequency) instead of just 
one (magnetic field strength) can be tuned. While optical 
FR have been employed to create molecules in cold gases 
via photoassociation (PA) ^3 and to tune the scattering 
length 0, they have not yet been used to coherently 
create molecules except for the recent work of Ref . Ql • 

In this Letter, we propose to employ optical FR to 
create weakly bound ground state molecules (in singlet 
and triplet ground state potentials both labelled 'ground 
state' in the following). In analogy to magnetic FR, we 
envisage a scheme of adiabatically ramping over the res- 
onance (cf. Fig. 0. The resulting wave function has 
components on both electronic ground and excited states 
with the latter being subject to spontaneous emission 
losses. In a second step, the laser field therefore needs 
to be switched off. This corresponds to projecting the 
wave function onto the field-free eigenstates. The goal is 
to sweep intensities and frequencies such that this pro- 
jection is predominantly onto the last bound level of the 
ground state, i.e. onto stable molecules. Our scheme is 



different from one-color PA [l^ which populates excited 
state levels. It is rather similar in spirit to Stimulated 
Raman Adiabatic Passage in that population of the 
excited state is minimized using a two-photon transition. 
It differs from two-color PA Q since the sudden switch-off 
breaks the symmetry of the coupling between the bound 
molecular level and the trap (or continuum) state. We 
show that for sufficiently tight traps, intensity and fre- 
quency of the field can be tuned such that spontaneous 
emission losses are avoided while adiabaticity is retained. 
Such ramps can be realized experimentally employing 
acousto-optical modulators or diode lasers. Sufficiently 
tight confinement can be reached in microscopic dipole 
traps or deep optical lattices 

Our calculations are performed for ®^Rb. The gen- 
erality of the scheme is emphasized by employing both 
singlet and triplet ground state potentials. We consider 
two *^Rb atoms which collide in an isotropic harmonic 
trap and interact with a continuous wave (CW) laser 
field. The center of mass motion is decoupled, and the 
dynamics in the internuclcar distance R is governed by 
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FIG. 1: (Color online) Our proposed scheme for the cre- 
ation of molecules: Step (1) is an adiabatic ramp of the 
laser intensity (and possibly frequency), while step (2) is 
a sudden switch-off of the laser (wave functions shown for 
utr = 50 kHz). 
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the Hamiltonian 
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Hg(e) = T + Vg(e)(R) + (-)Vtr(R) IS the single chan- 
nel Hamiltonian with T the kinetic energy operator and 
Vg(e)(R) the ground (excited) state interaction potential. 

Vtr(R) = ^muj'^j.R is the potential of the dipole trap, 
and r the decay rate modelling spontaneous emission, 
m denotes the reduced mass and cjtr the frequency of 
the trap (cjtr — 2Tr x i>tr)- The frequency of the laser, 
ujl = Wo — Ai is red-detuned by relative to the 
atomic resonance at loq. In Eq. we invoke the dipole 
and rotating wave approximations (RWA) . The Rabi fre- 
quency n is then given by SI = EoD{R) ■ e « EqD ■ e, 
where Eq is the amplitude of the laser field, D{R) the 
dipole moment and e the polarization vector of the laser 
field. -D(R) • e* is approximated by its asymptotic value 
deduced from standard long range calculations |l6|. In 
Eq. we neglect the hyperfine structure. This is jus- 
tified for sufficiently detuning the laser from the atomic 
resonances (about 4 cm~^ or 120 GHz, the largest en- 
ergy difference between hyperfine levels is 7 GHz between 
= 1 and F = 2 for S^Si/a). The potentials Vg(e)(R) 
have been obtained by matching the results of ab initio 
calculations ^3 to the long-range dispersion potentials 

Ksy(R) = {C3/R^+)Ce/R^ + Cs/R^. The coefficients 
for the 5S + 5S and 55 -I- 5P asymptote are respectively 
found in Ref. 18] and Ref. jlj,]. The repulsive barrier of 
the ground state potentials has been adjusted to give a 
triplet (singlet) scattering length of 100 ao (90 ao). The 
Hamiltonian, Eq. is represented on a grid, employing 
a mapping procedure [20j which reduces the number of 
required grid points by a factor of 5 to 30. 

We proceed in two steps. First, we diagonalize the 
Hamiltonian, Eq. QJ, and obtain the field dressed eigen- 
states and eigenenergies as a function of laser intensity 
and frequency. The term —ihT/2 causes the Hamil- 
tonian to be non-Hcrmitian with complex eigenvalues, 
r is assumed to be independent of R which is con- 
sistent with the approximation D{R) w D. There- 
fore the imaginary part of the eigenvalues becomes T/2 
times the projection of the eigenfunction onto the ex- 
cited state [23 ■ In a second step, we solve the time- 
dependent Schrodinger equation to illustrate the creation 
of molecules. T is then set equal to its asymptotic value, 
V2rat with Tat = ft/Tat and rat(5S + 5P3/2) = 26.24 ns, 
Tat(5S + 5Pi/2) = 27.70 ns. 

The following calculations are performed for transi- 
tions between the triplet ground state a^E+(5S-|-5S) and 
the 0~(5S-|-5P3/2) excited state. Fig.|2K shows the bind- 
ing energy of the last bound level below the (5S-I-5S) 
asymptote as a function of laser intensity and detuning. 
The range of detunings is chosen around 4 cm"-'^, large 
enough to avoid hyperfine coupling, and small enough 




FIG. 2: (Color online) (a) Binding energy of the last bound 
level below the (5S4-5S) asymptote as a function of laser in- 
tensity and detuning, (b) Binding energy of the last bound 
levels n as a function of intensity. For / =19 kW/cm^ (vertical 
line) there are Ntound levels below the (5S+5S) asymptote. 

such that the resonances occur with excited state levels 
v' which have a good Franck-Condon overlap with the 
last bound ground state level. Two resonances are found 
within this range (v' — 40 at 4.225 cm^^ and v' — 41 
at 3.98 cm~^). In Fig. IJb, the energies of the four last 
bound levels below the (5S-I-5S) asymptote are plotted 
vs. laser intensity for a specific detuning. Resonances at 
about 2.5 kW/cm^, 16.5 kW/cm^ and 36.5 kW/cm^ are 
observed. At each resonance, the number of bound states 
is increased by one. Usually, only the detuning is var- 
ied in optical FR. The increase in the number of bound 
states can then be understood as follows: in the RWA 
ground and excited state potential cross and the excited 
state asymptote is at HA l above the ground state disso- 
ciation limit. Decreasing the detuning therefore pushes 
one more excited state level below this dissociation limit. 
The same happens as intensity is increased. It corre- 
sponds to the light shifts displacing the resonance po- 
sitions with increasing intensity (cf. Fig. Et)- To fur- 
ther illustrate this "creation" of bound levels, Fig. 
(middle) shows the projection onto the ground state of 
one field dressed wave function, |(3|Vn=8i)P fo'^ differ- 
ent intensities, i.e. different (n counts all eigenstates). 
At / = 0, |</'^=8i) coincides with the lowest trap state 
(Fig. n right). As the intensity is increased, the wave 
function is deformed and pushed toward shorter inter- 
nuclear distances such that it eventually resembles the 
wave function of the last bound level (Fig. ^ left) [i^l . 
The first step in our scheme is therefore a slow ramp 
in intensity (and possibly frequency) such that the wave 
function adiabatically follows the field-dressed eigenfunc- 
tions, |</?^^gi). In a second step, the field should be 
suddenly switched off projecting the field-dressed onto 
the field-free eigenfunctions. The probability to form a 
ground state molecule is then given by the projection 
of the field-dressed eigenfunction onto the last bound 
ground state level, Pmoi = |(</'fastl'i'n=8i)P (^Iso lower 
bound levels can contribute to molecule formation, but 
this is much less likely). The lifetime of the field-dressed 
eigenfunction, = r^t / (V^Pcxc) , is determined by its 
excited state component, poxc = l(e|'/'n=8i) That is. 
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FIG. 3: (Color online) Probability of molecule formation Pmoi 
(a) and excited state component Poxc (b) as function of laser 
intensity and detuning (z^tr ~ 250 kHz). 
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TABLE I: Vibrational periods of the lowest trap state (calcu- 
lated from its eigenenergy) for the triplet and singlet ground 
state potentials. 



^'moi corresponds to a gain while Poxc might lead to a 
loss. Both are shown in Fig. |2| as a function of laser in- 
tensity and detuning. Close to resonance with an excited 
state level (Al = —4.225 cm^^), at moderate intensi- 
ties (5 kW/cm^ < I < W kW/cm^) the projection onto 
the last bound level is 50% and higher, while the excited 
state population does not exceed 0.01, i.e. the lifetime 
of the field-dressed eigenstate is > 2 /is. This lifetime 
defines an upper limit for the time window within which 
the ramp across the resonance should be completed. The 
lower limit. Tad, is due to the requirement of adiabaticity. 
It is determined by the vibrational period of the lowest 
trap state, Tvib, which depends on the trap frequency and 
the interaction potential (cf. Table Q}. We can now esti- 
mate the timescales for our scheme: Assuming the ramp 
should be performed in a time Tad = 5 x Tvib to be adia- 
batic, spontaneous emission losses should be minimal for 
a trap frequency of ftr > 250 kHz (Tad ~ 210 ns). For 
i^tr ~ 50 kHz, Tad (~ 1.4 /is) and Ty are on the same 
order of magnitude, and spontaneous emission losses will 
play a role. 

To verify our conclusions from the time-independent 
picture, we have explicitly studied the creation of 
molecules solving the time-dependent Schrodinger equa- 
tion, 

^h^^\^it)) = H{m{t)), (2) 

with a Chebychev propagator. The time-dependence in 
H(^) is due to the linear ramp in Q (i.e. Eq or v^) and 
llJl (i-c. Ai), respectively. Fig. 0] shows the projection 
of the wave function ^{R;t) onto the last bound level 
and onto the lowest trap levels of the field-free Hamil- 
tonian vs. time. Also plotted is the overall loss due to 
spontaneous emission ( we assume that any population 
undergoing spontaneous emission is lost from the coher- 
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FIG. 4: (Color online) (c)-l-(d); Projection of the time- 
dependent wave function onto the last bound ground state 
level (solid red line) and the first trap states T (dashed lines) 
of the bare Hamiltonian. Also shown is the total popula- 
tion |{*(t)|*(t))p (dotted line). (a) + (b): Variation of the 
laser field strength E{t) (solid line) and detuning Al (dashed 
line). The maximum field intensity is / = 8 kW/cm^ for 
v = 250 kHz and / = 5 kW/cm^ for // = 50 kHz. 



ence of the scheme, and possibly from the trap). We 
combined a ramp in intensity with a ramp in frequency, 
both ramps are performed within 85 ns and 540 ns, re- 
spectively, i.e. 2 X Tvib. This turned out to be sufficient to 
retain adiabaticity. While this result might be surprising 
at first glance, it reflects that the main source of nonadi- 
abaticity is Rabi cycling which in turn is suppressed by 
spontaneous emission. Simulations without spontaneous 
emission showed that each ramp time had to be at least 
5 X Tvib for the ramp to be adiabatic. Fig.^J: shows that 
for I'tr — 250 kHz, almost 50% of the population can be 
converted into ground state molecules, while about 10% 
are lost. The remaining population is distributed over the 
lowest trap states. For i'tr — 50 kHz (Fig. 011) the losses 
are somewhat higher at 24%, but the conversion proba- 
bility still reaches almost 30%. Conversion probabilities 
higher than 50% could be obtained for even tighter traps 
allowing for faster ramps. 

We found that combining a ramp in intensity with a 
ramp in frequency is the most efficient way to create 
molecules. The initial ramp in intensity is performed 
with the laser frequency tuned in between two reso- 
nances. Such a ramp deforms the wave function already 
considerably while keeping the excited state population 
(cf. Fig. Ob) and hence the spontaneous emission loss ex- 
tremely small. In a second step the frequency is ramped 
toward the resonance. This ensures a maximum overlap 
with the bound state wave function (cf. Fig.|3Ji). Ramp- 
ing up the intensity with the laser tuned in between two 
resonances is not advantageous to create molecules. The 
overlap with the trap levels exceeds in this case by far 
the overlap with the bound state, i.e. one mainly excites 
higher trap states. If the intensity is ramped up with 
the laser tuned close to resonance, due to the stronger 
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P with u = 250 kHz 
P with iy = 50 kHz 


0.475 
0.288 


0.269 
0.203 


0.389 
0.246 


0.466 
0.284 



TABLE H: Final probability of ground state molecule forma- 
tion, P = |('I'(tflnai)|<(5g''"'°>|^, for different switch-off times of 
the laser 
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coupling with the excited state, the ramp needs to be 
slower to be adiabatic, and spontaneous emission losses 
are larger (cf. Fig.OJ. 

Finally, since in the experiment the laser cannot be 
switched off instantly, we address the question of a fi- 
nite switch-off time after the ramp. (cf. Table ITl|l . 
As can be seen in Tabic HU the field should ideally be 
switched off within 1 ns. The probability of molecule 
formation is then reduced by less than 1%. Switching off 
the field within 10 ns causes about one third of the cre- 
ated molecules to be lost, i.e. to be transferred back into 
a pair of trapped atoms. Due to the shorter timescales, 
this problem is more severe in tighter traps. 

We have performed the same set of calculations for a 
transition between the Xii;+(5S+5S) and 0+(5S+5Pi/2) 
potentials, which differ markedly from the previous ones. 
The common point is the asymptotic behaviour in 

the excited 0+(5S+5Pi/2) potential, also providing ex- 
tremely long-range levels with a large Franck-Condon 
overlap with the lowest trap state(s). We therefore find 
exactly the same pattern of binding energies, probabil- 
ity of molecule formation and excited state population 
as shown in Figs. I2I4I 25j . 

To summarize we have shown that in tight traps, 
loosely bound ground state molecules can be created effi- 
ciently and without loss of coherence by adiabatic ramp- 
ing over an optical FR. Both detuning and intensity are 
varied within an asymmetric scheme which involves first 
adiabatic following and then a fast switch-off. Three- 
body effects are neglected in our model: it is therefore 
applicable to low densities or ideally to a Mott-Insulator 
state |15|. We assumed the trap to be isotropic. However 
we expect our scheme to work also for anisotropic traps 
as long as no new timescale is introduced, i.e. as long 
as the lowest trap frequency is > 50 kHz. We point out 
that in shallow traps the creation of molecules by opti- 
cal FR does not seem to be feasible. In that case non- 
adiabatic processes or adiabatic schemes employing short 
laser pulses jl^l should be considered. For tight traps 
by contrast, calculations performed for rubidium atoms 
demonstrated conversion probabilites up to 50%. The 
presented scheme should work for any optical transition 
which is also efficient in photoassociation, involving long- 
range wells or resonant coupling in the excited state [21| . 
The range of possible applications therefore extends well 
beyond homonuclear systems and alkali atoms. 

We would like to thank P. Naidon, P. Grangier, J. 
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